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Abstract 


We have studied the transverse-momentum (pj) dependence of the inclusive J/xff production in p-Pb 
collisions at = 5.02 TeV, in three center-of-mass rapidity (ycms) regions, down to zero pj. 
Results in the forward and backward rapidity ranges (2.03 < jcms < 3.53 and —4.46 < ycms < —2.96) 
are obtained by studying the i/xj/ decay to while the mid-rapidity region (—1.37 < ycms < 

0.43) is investigated by measuring the e+e^ decay channel. The pj dependence of the J/xj/ production 
cross section and nuclear modification factor are presented for each of the rapidity intervals, as well as 
the J/i// mean pj values. Forward and mid-rapidity results show a suppression of the J/y/ yield, with 
respect to pp collisions, which decreases with increasing pj. At backward rapidity no significant J/y/ 
suppression is observed. Theoretical models including a combination of cold nuclear matter effects 
such as shadowing and partonic energy loss, are in fair agreement with the data, except at forward 
rapidity and low transverse momentum. The implications of the p-Pb results for the evaluation of 
cold nuclear matter effects on J/y/ production in Pb-Pb collisions are also discussed. 
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The suppression of charmonia, bound states of c and c quarks, and in particular of the J/i/r state, has 
long been proposed as a signature for the formation of a plasma of quarks and gluons (QGP) Q in 
ultrarelativistic nucleus-nucleus collisions. However, it was soon realized that charmonium production 
can also be modified by nuclear effects not necessarily related to QGP formation Q. These so-called cold 
nuclear matter (CNM) effects can be investigated by studying charmonium production in proton-nucleus 
(p-A) collisions as confirmed by the analysis of results obtained by several fixed-target (SPS |[^Q, 
HERA Q and Tevatron Q) and collider (RHIC |[7| and EHC @1^) experiments. 


Theoretical models have studied the production of charmonium in p-A collisions and the effects of the 
surrounding cold nuclear medium by introducing various mechanisms which include nuclear shadowing, 
gluon saturation, energy loss and nuclear absorption. Models 110 -121 inspired by Quantum ChromoDy- 
namics (QCD) describe charmonium production as a two-step process, with the cc pair created in a hard 
parton scattering, followed by its evolution into a bound state with specific quantum numbers. The pair 
creation is sensitive to the Parton Distribution Eunctions (PDEs) in both colliding partners and, at high 
energy, occurs mainly via gluon fusion. Although PDEs are known to be modified in a nuclear environ¬ 
ment, information on the dependence of such modifications on the fraction x (Bjorken-x) of the nucleon 
momentum carried by the gluons and on the four-momentum squared transferred in the scattering 
is still limited 1 13-{^|. Charmonium production measurements can therefore provide insight into the 
so-called nuclear shadowing, i.e., on how the nucleon gluon PDEs are modified in a nucleus. 
Modifications of the initial state of the nucleus are also addressed by approaches assuming that at suffi¬ 
ciently high energies, when the quark pair is produced from a dense gluon system carrying small x-values 
in the nuclear target, a coherent effect known as gluon saturation sets in. Such an effect can be described 
by the Color Glass Condensate (CGC) effective theory, which is characterized by a saturation momentum 


scale (2i)- When combined with a specific quarkonium production model 116 


171, it is able to provide 


predictions for charmonium production in p-A collisions. In the context of shadowing and CGC models, 
a measurement of the charmonium yield as a function of transverse momentum (px) and rapidity (y) is 
important as it gives access to specific ranges of values of the gluon x and/or Q^. 

In addition to these purely initial state effects, both the incoming partons and the cc pair propagating 
through the nucleus may lose energy by gluon radiation at the various stages of the charmonium forma¬ 
tion process 118|. The interference of gluons radiated before and after the hard production vertex can 
lead to coherent energy loss effects, expected to induce a modification of the charmonium kinematic 
distributions | |T^ . 

Einally, while travelling through nuclear matter, the evolving cc pair or, if crossing times are sufficiently 
large, the fully formed resonance, may break-up into open charm meson pairs. Although this mecha¬ 
nism, known as nuclear absorption, plays an important role at lower collision energies Q, at the EHC 
the contribution of this effect to the production cross section is expected to be small, due to the very short 
crossing time of the pair through the nuclear environment. 


Understanding the role of the cold nuclear matter effects outlined above is essential to further our knowl¬ 
edge of various aspects of the physics of strong interactions, and it is crucial for the interpretation of the 
results on charmonium production in heavy-ion collisions, where the formation of a QGP is expected. 
In such a hot and dense deconfined medium the color screening mechanism (the QCD analogue of the 
Debye screening in QED) can prevent the formation of the heavy-quark bound states, leading to a sup¬ 
pression of quarkonium production |T|. In addition, at EHC energies, the large charm quark density 
may lead to a (re)generation of charmonium by (re)combination of charm quarks in the QGP 

phase and/or when the system cools down and the formation of hadrons occurs. This effect enhances 
charmonium production and is expected to be particularly sizeable at low px- In heavy-ion collisions, 
a superposition of hot and cold nuclear matter effects is expected, and a quantitative evaluation of the 
latter is an important prerequisite for a detailed understanding of the former. At lower energy, both 
at SPS 122-241 and RHIC [25 261, a suppression of J/i/r production, in addition to the CNM effects 
estimated from p-A(d-A) collisions, was indeed observed. 
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A suppression of J/i/r production has been measured in Pb-Pb collisions at the LHC [2^ 311. It was quan¬ 
tified via the nuclear modification factor, i.e., the ratio of the Pb-Pb yields with respect to those measured 
in pp at the same energy, scaled by the number of binary nucleon-nucleon collisions. The suppression 
has been found to be stronger at forward rapidity and at high px 130 311, in agreement with expectations 
from (re)combination models. Similar to the lower energy experiments, accurate measurements in p-A 
collisions are needed to quantitatively assess the contribution of hot and cold nuclear matter effects in 
Pb-Pb. 


The first measurements of inclusive J/i/r production in p-Pb collisions at the LHC at -^^nn = 5.02 
TeV 1^1^ have shown a sizeable suppression, with respect to binary-scaled pp collisions, at forward ra¬ 
pidity (p-going side) and no suppression at backward rapidity (Pb-going side). The nuclear modification 
factors are in fair agreement with models based on nuclear shadowing ||^[^. Calculations including 


a contribution from coherent energy loss 119| also reproduce the data. Corresponding measurements for 


the less strongly bound V^(25') charmonium state are presented in |34|. In addition, an extrapolation to 


Pb-Pb collisions of the J/i/r suppression measured in p-Pb showed that the effects observed in Pb-Pb 
cannot be ascribed only to CNM Q. 

In this situation, a study of the transverse-momentum dependence of J/i// production at LHC energies for 
various rapidity regions is particularly interesting in order to: (i) reach a deeper understanding and better 
quantify the complicated interplay of CNM effects, which are expected to exhibit a well-defined kine- 
matical dependence | ^35[36 1; (ii) determine if the differential features of the Pb-Pb results that suggest 
the presence of (re)combination effects are still present when the contribution of CNM is considered. 


In this paper, we present ALICE results on the transverse-momentum dependence of the inclusive J/yr 
production in p-Pb collisions at = 5.02 TeV, measured in three center-of-mass rapidity CVcms) 

ranges: backward (—4.46 < ycms < —2.96), mid- (—1.37 < ycms < 0.43) and forward (2.03 < ycms < 
3.53). The data are from the 2013 LHC p-Pb run. 

At mid-rapidity, J/i/r are reconstructed in the e+e^ decay channel with the ALICE central barrel detectors, 
covering the pseudorapidity range |T]iab| <0.9. Lor the backward and forward rapidity analysis, J/i/r are 
detected, through their decay channel in the muon spectrometer, in the pseudorapidity range 

-4 < T]iab < -2.5. 


Due to the energy asymmetry of the LHC beams (Lp = 4 TeV and Lpb = 1.58 • Apb TeV, where Apb= 
208 is the Pb atomic mass number), the nucleon-nucleon center-of-mass is shifted, with respect to the 
laboratory frame, by Ay = 0.465 in the direction of the proton beam. Since data were collected in 
two configurations, interchanging the direction of the proton and the Pb beams in the LHC, the muon 
spectrometer acceptance covers the forward and backward ycms regions quoted above, where positive 
(negative) rapidities refer to the direction of the proton (Pb) beam. In the following, the notation p-Pb 
(Pb-p) will refer to the first (second) configuration. 


Lor the dielectron analysis, the central barrel detectors used for the reconstruction are the Inner 
Tracking System (ITS) | [37] | and the Time Projection Chamber (TPC) |381. The ITS contains six cylin¬ 
drical layers of silicon detectors, with the innermost layer at a radius of 3.9 cm with respect to the beam 
axis and the outermost layer at 43 cm. This detector is used for reconstructing the primary interaction ver¬ 
tex as well as vertices from different interactions and secondary vertices from decays of heavy-flavored 
particles. The TPC has a cylindrical geometry with an active volume that extends from 85 to 247 cm 
in the radial direction and 500 cm longitudinally. It is the main central barrel tracking detector and also 
provides particle identification via the measurement of the specific energy loss (dL/dr) in the detector 
gas. 


The muon spectrometer |39| is the main detector used in the dim non analysis. It consists of a 3 T-m 
dipole magnet, coupled with a tracking and a triggering system. Between the interaction point and the 
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muon spectrometer, a ten interaction-length (Ai) front absorber filters out the hadrons produced in the 
interaction. Muon tracking is performed by means of five tracking stations, each one made of two planes 
of Cathode Pad Chambers. A 7.2 Aj iron wall, which stops secondary hadrons escaping the front absorber 
and low momentum muons, is placed after the tracking stations. It is followed by a muon trigger system, 
based on two stations equipped with Resistive Plate Chambers. A conical absorber made of tungsten, 
lead and steel protects the spectrometer against secondary particles produced by the interaction of large-T] 
primary particles in the beam pipe. In the dimuon analysis, the determination of the interaction vertex is 
provided by the two innermost Si-pixel layers of the ITS (Silicon Pixel Detector, SPD). 


For both analyses, timing information from the Zero Degree Calorimeters |40|, placed symmetrically 
at 112.5 m with respect to the interaction point, is used to remove de-bunched proton-lead collisions. 
Furthermore, two scintillator hodoscopes (VZERO) |411, with pseudorapidity coverage 2.8 < T]iab <5.1 
and —3.7 < T]iab < —1-7, are used to remove beam-induced background. More details on the AFICE 
apparatus can be found in 1391. 


A coincidence of signals in the two VZERO detectors provides the minimum bias (MB) trigger, which 


has a > 99% efficiency for selecting non single-diffractive p-Pb collisions |42|. While the dielectron 
analysis is based on MB-triggered events, the study of J/i// in the ^ decay channel relies on a dimuon 
trigger which requires, in addition to the MB condition, the detection of two opposite-sign tracks in the 
trigger system. The dimuon trigger selects two muon candidates with transverse momenta larger 
than 0.5 GeV/c. The trigger threshold is not sharp, and the single muon trigger efficiency reaches its 
plateau value (~ 96%) at px./j ~ 1-5 GeV/c. The dielectron analysis was performed on a data sample 
corresponding to the p-Pb configuration, with an integrated luminosity = 51.4 ± 1.9 while 

for the dim lion analysis the corresponding values are 5.01 ±0.19 nb^^ for p-Pb and 5.81 ±0.20 nb^^ 


for Pb-p (the quoted uncertainties are systematic) |431. 


The dielectron analysis is based on 1.07x10^ events, collected with a low MB interaction rate (~10 
kHz), with a negligible amount of events having more than one interaction per bunch crossing (pile- 
up events). The interaction vertex is required to lie within ±10 cm from the nominal collision point 
along the beam axis, in order to obtain a uniform acceptance of the central barrel detector system in 
the fiducial range Ifiiabl < 0.9. Electron candidates are selected with criteria very similar to those used 
in previous analyses of pp collisions at ±5 =7 TeV [441 and Pb-Pb collisions at i/^nn = 2.76 TeV 
pO) . To ensure a uniform tracking efficiency and particle identification resolution in the TPC, only 
tracks within |T]iab| < 0.9 are used. Electron identification is performed using the TPC, as shown in 
Fig. by requiring the dE/dr signal to be compatible with the electron assumption within 3a, where 
a denotes the resolution of the AE/Ax measurement. Furthermore, the TPC tracks that are compatible 
with the pion and proton assumptions within 3.5a are rejected. A slightly looser rejection condition 
(3a) is applied when considering tracks corresponding to dielectron candidates with > 5 GeV/c 
in order to enhance the statistics. A cut on the transverse momentum (px e > 1.0 GeV/c) is applied 
to remove combinatorial background from low-momentum electrons. The efficiency loss induced by 
this cut amounts to only ~20%, due to the relatively large momentum of the J/i/r decay products. The 
electron candidates must have at least one hit in the innermost two layers of the ITS, thus rejecting a 
large fraction of background electrons from photon conversions. For dielectrons with px < 3 GeV/c the 
electron candidates are required to have a hit in the first layer, to further reduce background. The tracks 
are required to have at least 70 out of a maximum of 159 clusters in the TPC and a normalized to the 
number of clusters attached to the track smaller than 4. 


The J/i/r yields are obtained by counting the number of entries in the invariant mass range 2.92 < mg+e- < 
3.16 GeV/c^ after background subtraction. The J/i// radiative decay channel and the energy loss of the 
electrons due to bremsstrahlung in the detector material produce a long tail towards low invariant masses. 
A fit using a Crystal Ball (CB) | |45| function for the signal gives compatible values in Monte-Carlo 
(MC) and data (~20 MeV/c^ for the width of the Gaussian component of the CB). Taking into account 
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Fig. 1: Charged particle specific energy loss (dL/dr) as a function of momentum, as measured in the TPC in p-Pb 
collisions. The black lines are the corresponding Bethe-Bloch parametrizations for the various particle species. 



(GeV/c2) 


Fig. 2: Opposite-sign dielectron invariant mass spectra (blue symbols) for various pj intervals, compared to the 
background (black curve) estimated through mixed events. The background is scaled to match the data in the mass 
ranges 2.0 < m^+^- < 2.5 GeV/c^ and 3.2 < m^+g- < 3.7 GeV/c^. 

such a mass resolution and the presence of the bremsstrahlung tail, 67 — 73% of the signal, depending 
on pt> falls within the counting window. The background shape is obtained from event mixing. Event 
mixing is performed by pairing leptons from different events having similar global characteristics such 
as the primary-vertex position and the track multiplicity (the result being quite insensitive to the rapidity 
range, either forward or central, chosen for the multiplicity measurements). The mixed-event background 
is then scaled to match the same-event opposite-sign distribution in the mass ranges 2.0 < mg+e- < 
2.5 GeV/c^ and 3.2 < nje+e- <3.7 GeV/c^ (the contribution of the bremsstrahlung tail in the former 
range and of the '/^(2S) in the latter are negligible). Consistent results are found when the same-event 
like-sign distributions are used, instead of event mixing, to estimate the background. The systematic 
uncertainty on the signal extraction comes from the variation of the mass range where the normalization 
of the mixed-event background shape is performed and from the choice of the mass window where the 


5 



























Rapidity and transverse-momentum dependence of the inclusive J/i/r RpPb 


ALICE Collaboration 


signal is counted. The signal extraction has been performed in five transverse-momentum bins, pj <13, 
1.3 < pt < 3, 3 < pt < 5, 5 < pt < 7 and 7 < px < 10 GeV/c. The J/ip counts in these bins vary 
from 25 to 132, with a significance, computed in the 2.92 < < 3.16 GeV/c^ mass region, ranging 

from 4.6 to 8.7. An analysis of the px-integrated data sample, using the procedure detailed above, gives 
465 ±37(stat.) ± 16(syst.) J/ip signal counts. The systematic uncertainty on the signal extraction is 
largest at low px (10% for px < 1.3 GeV/c and 12% for 1.3 < px < 3 GeV/c), due to a less favorable 
signal over background ratio, and decreases to ~5.5-8.4% in the other three px bins. Figurej^shows the 
invariant mass distributions for the opposite-sign dielectrons compared with the mixed-event background 
for the different intervals of px. 


The dimuon analysis is performed as detailed in Q, and is shortly summarized hereafter. Data were 
collected with the dimuon trigger, and the MB interaction rate (up to 200 kHz) was much higher than 
in the sample used for the dielectron analysis. This leads to a ~2% interaction pile-up probability. 
However, the probability of having more than one dimuon in the same bunch crossing satisfying the 
trigger condition is negligible. Muon candidate tracks are reconstructed in the tracking system by using 
the standard reconstruction algorithm 1441. The quality of the tracks is ensured by requiring the single 
muon pseudorapidity to be in the range —4 < Tjiab./i < —2.5, in order to remove particles at the edges 
of the muon spectrometer acceptance. In addition, a cut on the radial coordinate of the track at the 
end of the front absorber (17.6 < Rabs < 89.5 cm) is performed, ensuring rejection of muons crossing 
its high-density part, where energy loss and multiple scattering effects are more important. The tracks 
reconstructed in the tracking system that are not matched to a corresponding track in the triggering 
system are rejected | |44| . Finally, the reconstructed dimuons are required to be in 2.03 < ycms <3.53 
(—4.46 < ycms < —2.96) for the forward (backward) rapidity analysis. The number of J/ip is extracted 
in transverse-momentum bins, in the range px < 15 GeV/c, through fits to the invariant mass spectra of 
opposite-sign dimuons. The spectra are fitted with a superposition of background and resonance shapes. 
The background is described with a Gaussian function with a mass-dependent width or, alternatively, 
with an exponential function times a fourth-order polynomial function. For the J/ip shape an extended 
Crystal Ball function, which accommodates a non-Gaussian tail both on the right and on the left side of 
the resonance peak, is adopted. Alternatively, a pseudo-Gaussian function | [46| is used, corresponding to 
a Gaussian core around the J/ip pole, and tails on the right and left side of it, parameterized by varying 
the width of the Gaussian as a function of the mass. The value of the J/ip mass and its width (a) at the 
pole position are free parameters of the fit. The mass coincides with the PDG value within less than 5 
MeV/c^ and the width is ~70 MeV/c^, slightly increasing with px, due to a small relative decrease in the 
tracking resolution for harder muons. Although the signal over background ratios, calculated for a ±3a 
interval around the resonance peak, are relatively large (ranging from 1.4 to ~ 6 moving from low to high 
Px), the parameters of the tails of the J/ip distributions cannot be reliably tuned on the data (in particular 
at large px, where statistics is limited), but are fixed, for each px bin, to the values extracted from fits 
to reconstructed samples from a signal-only MC generation. The contribution of the V^(2S) resonance 
is also included in the fitting procedure, even if its influence on the determination of the J/ip yield is 
negligible. Finally, all the fits are performed in two different invariant mass ranges, either 2 < < 5 

GeV/c^ or 2.2 < mfxfi < 4.5 GeV/c^. Examples of fits to the invariant mass spectra, in the px bins under 
study, are shown in Fig.[^ 


For each px bin, the number of J/ip is evaluated as the average of the integrals of the resonance functions 
obtained in the various fits. The RMS of the corresponding yield distributions (0.2 — 3%, depending on 
Px) provides the systematic uncertainty on the signal extraction. Additional sets of tails, obtained from 
the MC, but referring to other ycms and px phase space regions, have also been tested and the dependence 
of the extracted yields on the variation of the tails (2%) is included in the systematic uncertainty on 
the signal extraction. As a function of px, the number of J/ip in the p-Pb (Pb-p) configuration ranges 
between ~16100 (~16000) in the most populated bin (1 < px < 2 GeV/c) and less than ~900 (~300) 
in the highest px bin (10 < px < 15 GeV/c). 


6 




Rapidity and transverse-momentum dependence of the inclusive J/i/r RpPb 


ALICE Collaboration 



Fig. 3: The opposite-sign dimuon invariant mass spectra for the various px bins, relative to the p-Pb data sample 
(blue symbols). The fits shown in this Figure (blue curves) were performed by using the sum of extended Crystal 
Ball functions for the J/i// and V^(2S) signals, and a variable width Gaussian for the background. The signal and 
background components are shown separately as red curves. 


The J/i/r yields are then corrected for the product of acceptance times efficiency (Axe), evaluated by 
means of a MC simulation. J/y/ production is assumed to be unpolarized, as motivated by the small 
degree of polarization measured in pp collisions at = 7 TeV 147-49|. In the e+e^ decay channel, 
A X £ is calculated using a MC simulation where J/i/r are injected into p-Pb collisions simulated with 
HIJING |50|. The decay products of the J/i/r are then propagated through a realistic description of 
the ALICE set-up, based on GEANT3.21 | [5T[ , taking into account the time evolution of the detector 
performance. Einally, J/y/ candidates are reconstructed with the same procedure applied to data. The pj- 
integrated A x e factor amounts to 8.9%. Its px-dependence exhibits a minimum (~7.5%) around px = 
2 GeV/c, due to the kinematical acceptance, and it reaches ~12% at high px. The integrated value of A x 
£ is affected by a 3% systematic uncertainty related to the choice of the J/y/ pj- and y-distributions used 
in the MC simulation. This value is obtained using as input several distributions, determined by varying 
within uncertainties the differential spectra extracted from the ALICE p-Pb data themselves. Eor px- 
differential studies, the values of A x £ are found to be sensitive only at a sub-percent level to the adopted 
input Px- and y-distributions. A further small systematic uncertainty reaching 1.5% in the highest px 
interval and related to the statistical uncertainty of the MC sample is also introduced. The systematic 
uncertainty on the dielectron reconstruction efficiency is strongly dominated by the particle identification 
uncertainty and amounts to 4%. It was obtained by comparing the single track reconstruction efficiency 
for topologically identified positrons and electrons from photon conversions with the corresponding MC 
quantities. In the dimuon analysis, the J/y/ A x £ is obtained with a MC simulation, by generating 
signal-only samples, tracking them in the experimental set-up modeled with GEANT3.21 and using 
the same reconstruction procedure applied to data. The use of a pure signal MC is justified, since the 
tracking efficiency does not show a dependence on the hadronic multiplicity of the collision. A realistic 
description of the set-up is adopted, including the time evolution of the efficiencies of tracking and 
triggering detectors. As for the dielectron analysis, the differential distributions used as an input to the 
MC are tuned directly on the data. The J/ip A x £ values, integrated over px, are 25.4% and 17.1% 
for p-Pb and Pb-p respectively Q, and exhibit a dependence on transverse momentum, being of the 
order of ~24% (~16%) for p-Pb (Pb-p) at low px and reaches ~50% (~35%) in the highest px bin 
(10 < Px < 15 GeV/c). The systematically lower A x £ values in Pb-p reflect the smaller detector 
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Source 

*^pPb ’ RpPb 

-1.37<ycms <0.43 

^pPb ’ RpPb 
2.03<ycms <3.53 

^Pbp ’ RPbp 
-4.46< ycms <-2.96 

Uncorrelated 

Tracking efficiency (p+p^) 

_ 

4 

6 

Trigger efficiency (p+p^) 

- 

2.7 -4.1 

2.7 -4.1 

Matching efficiency (p+p^) 

- 

1 

1 

Reconstruction efficiency (e+e^) 

4 

- 

- 

Signal extraction 

5.5 - 12.6 

2 - 2.5 

2 - 3.6 

MC input 

0.3 - 1.5 

0.1 -0.4 

0.1 - 1.4 

Opp 

4.8 - 15.7 

5.2 - 9.2 

5.2 - 9.2 

Partially correlated 
app"^ (corr. vs y and px) 


2.8 - 5.9 

2-5.6 

Correlated 

B.R. (J/ip^Z+Z^) 

1 

1 

1 

^inticorr. vs. px, uncorr. vs. y) 

3.3 

3.4 

3.1 

i:^„,(corr. vs. y and px) 

Opp 

1.6 

16.6 

1.6 

5.2 

1.6 

5.2 


Table 1: Systematic uncertainties (in percent) on the measurement of inclusive J/y/ cross sections and nuclear 
modification factors. For px-dependent uncertainties, the minimum and maximum values are given. The degree of 
correlation (uncon'elated, partially correlated, correlated) refers to the px-dependence, unless specified otherwise. 
It cannot be excluded that a degree of correlation, difficult to quantify, is present also in uncertainties cun'ently 
labelled as uncorrelated. Uncertainties on ££int and branching ratios are relevant for cross sections, while those on 
app"^ contribute only to the uncertainty on the nuclear modification factors. uncertainties are split into two 
components, respectively uncorrelated and correlated between p-Pb and Pb-p, as detailed in | |43) . 

efficiency in the corresponding data taking period. The systematic uncertainty on the integrated Axe due 
to the input shapes is 1.5% for both p-Pb and Pb-p, and has been estimated using various distributions 
obtained from data and corresponding to smaller intervals in y, px and centrality (see |[8| for details). 
For px-differential studies, the corresponding uncertainties are below 1.5%. The uncertainty on the 
dimuon tracking efficiency amounts to 4% (6%) for p-Pb (Pb-p) and is taken as constant for the full pj 
range. It is evaluated by combining the uncertainties on single muon tracking efficiencies, considered 
as uncorrelated. The efficiency of each tracking plane is obtained using the redundancy of the tracking 
system (two independent planes per station) and then single muon efficiencies for the full tracking system 
are calculated according to the tracking algorithm | [52| . Their uncertainty is determined by comparing the 
efficiency obtained with tracks from MC and real data. The systematic uncertainty on the dimuon trigger 
efficiency includes: (i) a contribution due to the uncertainty in the evaluation of the trigger detector 
efficiency (~ 2%, independent of px); (h) a 0.5 — 3% px-dependent contribution (2% for the integrated 
efficiency), related to small differences in the trigger response function between data and MC in the 
region close to the trigger threshold; (iii) a 0.5 — 3.5% px-dependent contribution due to a small fraction 
of opposite-sign pairs which were misidentified as like-sign by the trigger system. Finally, a ~1% 
uncertainty, independent of px, is included, due to the choice of the value of the cut applied to the 
matching of tracks reconstructed in the muon tracking and triggering systems. 

The differential cross section for inclusive J/i/r production is defined as: 


__ (Vj/y,(Ay,Apx) _ 

dydpx . (A X £)(ax,Apt) ' B-R- {3/W ^ /+/-)• Ay • Apx 


( 1 ) 


where Aj/^(Ay, Apx) is the number of J/y/ for a given Ay and Apx interval. The branching ratio to dilep 
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tons, B.R.(J/i/r —)■ /+Z ), is 5.94 ±0.06% (5.93 ±0.06%) for the dielectron (dimuon) decay |531. The 
integrated luminosity, 


, is the ratio between A^MB^ the number of MB collisions, and CjM®, the corre- 


spending cross section, measured in a van der Meer scan to be 2.09 ± 0.07 b for the p-Pb configuration 
and 2.12 ± 0.07 b for the Pb-p case |43|. The luminosity is also independently determined by means 
of a second signal based on a Cherenkov counter |39|, as described in |43|. The two measurements 
differ by at most 1 % throughout the whole data-taking period and such a value is quadratically added 
to the luminosity uncertainty. Finally, since the dimuon analysis is based on a sample of Adimu dimuon 
triggered events, the number of equivalent MB collisions is computed as Amb = E • ADIMU^ where F is a 
factor accounting for the probability of having a dimuon trigger when the MB condition is satisfied and 
for the small (~ 2%) pile-up probability in the corresponding data sample. The systematic uncertainty 
on this quantity, quadratically added to the other luminosity uncertainties, is 1% and originates from the 
comparison between the different approaches used for its evaluation @ . A summary of the systematic 
uncertainties can be found in Table [T] The differential inclusive J/i/r cross sections are shown in Fig. 
in the ranges pj < 10 GeV/c for the dielectron analysis and px < 15 GeV/c for the dimuon analysis. The 
numerical values can be found in Table lU 


For the dielectron analysis, the px-integrated cross section was also determined, obtaining 
dapp^/dy(-1.37 < ycms <0.43) = 909 ±78(stat.) ±71(syst.)pb. 

The corresponding px-integrated cross sections for the dimuon analysis were published in ||8|. 
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Fig. 4: px-differential inclusive J/y/ cross sections for the various rapidity regions under study. The vertical error 
bars correspond to the statistical uncertainties, while open boxes represent the uncorrelated systematic uncertainties 
and the shaded boxes the quadratic sum of the fully and partially correlated ones. The numerical values can be 
read in Table|^ The horizontal bars correspond to the widths of the px bins. 


Starting from the px-differential J/ip cross sections it is possible to evaluate, as additional information, 
the mean px ((px)) for the various y-ranges, by means of fits based on the empirical function: 


o r>ppt, 

dydpx 
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PT 

(GeV/c) 

d^OppJf/dydpT 

(Alb/(GeV/c)) 

PT 

(GeV/c) 

^pPb 

d^CJpp'^/dydpx (interpol.) 
(Alb/(GeV/c)) 

4.46 < Jems 2.96 fx ) 

[0;i] 

97.7±2.0±7.2±3.5 

[0;i] 

0.96±0.02±0.09±0.03±0.06 

0.490±0.029±0.017±0.026 

[i;2] 

196.8±2.7±14.3±7.1 

[1;2] 

1.06±0.01±0.10±0.04±0.07 

0.892±0.048±0.030±0.046 

[2; 3] 

159.6±2.1±11.6±5.8 

[2; 3] 

1.11 ±0.01 ±0.10±0.04±0.07 

0.693±0.036±0.025±0.036 

[3;4] 

93.3±1.6±6.7±3.4 

[3;4] 

1.16±0.02±0.10±0.04±0.07 

0.388±0.021±0.012±0.020 

[4; 5] 

45.7±1.0±3.2±1.7 

[4; 5] 

1.17±0.02±0.11 ±0.03±0.07 

0.187±0.011±0.004±0.010 

[5; 6] 

22.1±0.5±1.6±0.8 

[5; 6] 

1.13±0.03±0.12±0.02±0.07 

0.094±0.007±0.002±0.005 

[6; 7] 

11.2±0.4±0.8±0.4 

[6; 8] 

1.27±0.03±0.14±0.08±0.08 

0.032±0.003±0.002±0.002 

[7; 8] 

5.7±0.3±0.4±0.2 




[8; 10] 

2.3±0.1±0.2±0.1 




[10; 15] 

0.33±0.03±0.03±0.01 




-1.37 <ycms < 0.43 (e+e ) 

[0;1.3] 

158±33±17±6 

[0;1.3] 

0.81±0.17±0.10±0.14 

0.94±0.07±0.16 

[1.3;3] 

211±33±26±8 

[1.3;3] 

0.64±0.10±0.09±0.11 

1.60±0.08±0.26 

[3; 5] 

126±15±9±5 

[3; 5] 

0.77±0.09±0.07±0.13 

0.79±0.05±0.13 

[5; 7] 

43.4±6.5±3.4±1.7 

[5; 7] 

0.89±0.13±0.13±0.15 

0.23±0.03±0.04 

[7; 10] 

10.2±2.4±1.0±0.4 

[7; 10] 

0.89±0.21±0.16±0.15 

0.06±0.01±0.01 

2.03 <ycms< 3.53 (p+p ) 

[0;i] 

78.8±1.5±4.6±3.1 

[0;i] 

0.61±0.01±0.05±0.02±0.04 

0.624±0.036±0.025±0.032 

[i;2] 

158.4±2.2±9.0±6.2 

[i;2] 

0.64±0.01±0.05±0.02±0.04 

1.197±0.064±0.046±0.062 

[2; 3] 

138.2±1.9±7.9±5.4 

[2; 3] 

0.68±0.01±0.05±0.03±0.04 

0.980±0.051±0.039±0.051 

[3;4] 

91.3±1.4±5.0±3.6 

[3;4] 

0.76±0.01±0.06±0.03±0.05 

0.579±0.032±0.022±0.030 

[4; 5] 

53.0±0.9±2.8±2.1 

[4; 5] 

0.87±0.02±0.07±0.02±0.06 

0.294±0.017±0.008±0.015 

[5; 6] 

29.7±0.6±1.6±1.1 

[5; 6] 

0.91±0.02±0.08±0.03±0.06 

0.156±0.011±0.005±0.008 

[6; 7] 

14.9±0.4±0.8±0.6 

[6; 8] 

0.98±0.02±0.09±0.05±0.06 

0.057±0.005±0.003±0.003 

[7; 8] 

8.3±0.3±0.5±0.3 




[8; 10] 

3.7±0.1±0.2±0.1 




[10; 15] 

0.77±0.03±0.05±0.03 





Table 2: Summary of the results on the inclusive J/y/ differential cross sections and nuclear modification factors 
for p-Pb collisions. The results of the cross section interpolation for pp collisions are also shown. For p-Pb cross 
section results, the first quoted uncertainty is statistical. The following uncertainties are systematic, the second 
one being pj-uncorrelated and the third one pj-correlated. For /?ppb the first quoted uncertainty is statistical. The 
following uncertainties are systematic, the second one being pp-uncorrelated. For dielectron results the third un¬ 
certainty is pT-correlated, while for dimuon results the third uncertainty is partially pp-correlated and the fourth is 
px-correlated. For the results on the interpolated pp cross section, the first quoted uncertainty combines statistical 
and pT-uncorrelated systematic uncertainties. For dielectron results the second uncertainty is px-correlated sys¬ 
tematic, while for dimuon results the second uncertainty is partially px-correlated, and the third is px-correlated. 
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where C, po and n are free parameters. The quality of the fits is satisfactory (;^^/ndf ~ 1) and the resulting 
(pt) values, computed for the measured pT ranges, are 

(pt)(— 4.46 < ycms < —2.96) = 2.47 ±0.01 (stat.) ±0.03(syst.)GeV/c 
(pt)(— 1.37 < ycms < 0.43) = 2.86±0.15(stat.) ±0.10(syst.)GeV/c 
(pt) (2.03 < ycms < 3.53) = 2.77±0.01(stat.)±0.03(syst.)GeV/c 


The quoted uncertainties were obtained by performing fits including only statistical (or uncorrelated 
systematic) uncertainties on differential cross sections. 

In order to perform a meaningful comparison of (px) results in the dielectron and dimuon analysis, the 
values from the dimuon analysis have also been extracted, with the same procedure detailed above, in 
the range px < 10 GeV/c, obtaining results which are smaller by less than 2% with respect to the full px 
range. It is found that (px) is larger at central rapidity. Furthermore, the (px) measured at forward ycms 
is significantly larger than at backward ycms- This difference, which could be partly due to the slightly 
different |y|-coverage, persists when (px) is calculated in the |ycms| region common to p-Pb and Pb-p 
(2.96 < |ycms| < 3.53). The values obtained in this case are 2.58 ±0.02(stat.) ±0.04(syst.)GeV/c and 
2.69 ±0.02(stat.) ±0.03(syst.)GeV/c, respectively at backward and forward ycms> and differ by ~ 2a. 

The J/ip nuclear modification factor RpPb is obtained as the ratio of the differential cross sections between 
proton-nucleus and proton-proton collisions, normalized to Apbi 

Apb • d2ap^'''/dydpx 



Since no pp data are available at ^/s = 5.02 TeV, the d^app'^/dydpx reference cross sections were ob¬ 
tained by means of an interpolation/extrapolation procedure. For the dielectron analysis, the starting 
point of the interpolation procedure is the determination of da/dy for inclusive J/ip in pp collisions at 
ycms ~ 0 and = 5.02 TeV, carried out as for the analysis described in 130|. Available mid-rapidity data 
at ^/s = 0.2 154|, 1.96 |(^, 2.76 |56| and 7 TeV |44| are interpolated using several empirical functions 


(exponential, logarithmic and power-law, covering in this way the various possibilities for the curvature 
of the \/i-dependence) obtaining da/dy = 6.19 ± 1.03 /ib. Even if the ycms range covered in this anal¬ 
ysis is shifted by 0.465 units with respect to mid-rapidity, the rapidity-dependence of the cross section 
is negligible compared to the uncertainty on the interpolation procedure. Then, a method similar to the 


one in |57| is applied to derive the px-differential cross section. It is based on the empirical observa¬ 
tion that pp and pp results on differential spectra obtained at various collision energies and in different 
rapidity ranges [44 48 54 55 ^ exhibit scaling properties when plotted as a function of px/(7’t)- The 


normalized spectra, with the statistical and the bin-by-bin uncorrelated systematic uncertainties added 
in quadrature, can be fitted with a one-parameter function described in | [57| . The px-differential cross 
sections at mid-rapidity and = 5.02 TeV can then be obtained by rescaling the fitted universal dis¬ 
tribution using the previously estimated da/dy and its corresponding (px)- The latter value is obtained 
by an interpolation of the energy-dependence of (px) values evaluated fitting the available experimental 
mid-rapidity results [44 54 ^ with exponential, logarithmic and power-law functions. One obtains in 


this way, in the range px < 10 GeV/c, (px) = 2.81 ±0.10 GeV/c as an average of the results calculated 
with the various empirical functions. As outlined above for da/dy, the 0.465 y-unit shift of the data with 
respect to mid-rapidity has a negligible effect also on (px). 

For the dimuon analysis, thanks to the smaller uncertainties with respect to mid-rapidity results, an 
approach equivalent to that described in |59|, exclusively based on the AFICE data collected at ^/s = 
2.76 TeV and 7 TeV |[60| in 2.5 < ycms < 4, px < 8 GeV/c has been used. The reference cross 
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sections are obtained with a two-step procedure, corresponding to an energy interpolation followed by 
a rapidity extrapolation. In the first step, for each px bin, the d^app'^/dydpx values at = 2.76 and 
7 TeV are interpolated, using three different empirical functions (linear, power-law and exponential) to 
estimate the cross section values at ^ = 5.02 TeV. The central values are calculated as the average of the 
results obtained with the three functions, while the associated uncertainties come from the experimental 
uncertainties on the points used for the interpolation, added in quadrature to a contribution chosen as the 
maximum spread of the results from the different interpolating functions. In the second step, this result 
is extrapolated from 2.5 < ycms < 4 to the p-Pb and Pb-p ycms ranges, using the scaling factors for the 
px-integrated cross sections computed in |59|. Finally, since the LHCb Collaboration has shown that the 
J/ip px distributions slightly depend on ycms |481 in the rapidity range covered in the dimuon analysis, a 
px-dependent correction tuned on these data (10% maximum at large px) is applied. 


The inclusive J/ip nuclear modification factor is shown in Fig.j^for the three rapidity regions under study. 
The numerical values of Rppb, as well as the results of the interpolation procedure for the estimate of the 
pp cross sections, can be found in Table]^ For the dimuon analysis, the evaluation of RpPb is restricted to 
Px < 8 GeV/c, the region covered by the pp measurements used in the evaluation of the reference cross 
sections. The sources of systematic uncertainties on RpPb and their values are summarized in Table[T] The 
terms related to the pp reference cross sections contribute to uncorrelated, partially or fully correlated 
uncertainties on Rppb, depending on their origin. In particular, for the dimuon analysis: (i) the statisti¬ 
cal and px-uncorrelated systematic uncertainties on the ^/s = 2.76 and 7 TeV pp data contribute to the 
uncorrelated uncertainty; (ii) the spread of the results obtained with various interpolating/extrapolating 
functions in ^ and ycms contribute to the partially correlated uncertainty; (iii) the yA'-correlated uncer¬ 
tainties between the ^/s = 2.76 and 7 TeV pp data contribute to the correlated uncertainty. At forward 
and mid-rapidity the J/ip RpPb shows a clear suppression at low px, vanishing at high px. At backward 
rapidity no suppression is present, within uncertainties. 


For the dielectron analysis, the px-integrated nuclear modification factor was also calculated, carrying 
out the signal extraction procedure on the px-integrated invariant mass spectrum. The obtained value 


^pPb = 0.71 ±0.06(stat.) ±0.13(syst.) 


is consistent with the forward rapidity (2.03 < ycms < 3.53) dimuon result, and smaller than the backward 
one (—4.46 < ycms < —2.96) by ~ 2a 

In Fig. 1^ predictions from various models are compared to the data. A calculation based on the next- 
to-leading order (NFO) Color Evaporation Model (CEM) for the prompt J/ip production and the EPS09 
shadowing parametrization | [^ reproduces within uncertainties the px-dependence and the amplitude 
of the suppression for px > 1.5 GeV/c in the three rapidity regions under study. The theoretical uncer¬ 
tainties arise from the uncertainties on EPS09 as well as on the values of charm quark mass and of the 
renormalization and factorization scales used for the cross section calculation. Data are also compared 
to two calculations based on a parametrization of experimental results on prompt J/ip production in pp 


collisions and including the effects of coherent energy loss 1351 in the cold nuclear medium. One of 
the calculations includes only coherent energy loss, while the other combines coherent energy loss with 
EPS09 shadowing. The uncertainty bands include, for the coherent energy loss mechanism, a variation 
of both the qo parameter (gluon transport coefficient evaluated at x = 0.01) and the parametrization of 
the production cross section. At forward rapidity the pure energy loss scenario predicts a much steeper 
px-dependence, while better agreement is found when the EPS09 contribution is included. However, at 
low Px, a discrepancy between data and both calculations is observed. Also at mid-rapidity the coherent 
energy loss model including the EPS09 contribution better describes the data, although the larger un¬ 
certainties prevent a firm conclusion. The same features can be observed at backward rapidity, where 
the calculation including coherent energy loss and shadowing agrees with the data in showing weak nu¬ 
clear effects on J/ip production. Finally, the results at central and forward rapidities are compared with 
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Fig. 5 : The J/y/ nuclear modification factor as a function of px at backward (top), mid (center) and forward 
(bottom) rapidities. Statistical uncertainties are represented by vertical error bars, while open boxes correspond 
to uncorrelated uncertainties and the shaded areas to uncertainties partially correlated in pj. The boxes around 
f^pPb = 1 show the size of the correlated uncertainties. The horizontal bars correspond to the widths of the px bins. 
Results from various models are also shown, including a pure shadowing calculation based on the EPS09 
parameterization, a CGC-inspired model and the results of the coherent energy loss calculation p5j , with or 
without the inclusion of an EPS09 shadowing contribution. 
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a prediction based on the CGC framework and using CEM for the prompt J/i/r production |36|. In the 
backward rapidity region, higher gluon v in the nucleus are probed and the CGC model is out of its range 
of applicability. The quoted uncertainties are related to the choices of and of the charm quark mass. 
While the model is in fair agreement with mid-rapidity data, it clearly underpredicts the J/y/ Rppt in the 
full pj range at forward rapidity. 


The theoretical calculations discussed above are carried out for prompt J/y/ (i.e., direct J/y/ and the 
contribution from Xc and yz/lS) decays), while the measurements are for inclusive J/y/ which include 
a non-prompt contribution from B-hadron decays. The contribution of the latter source to can 
be evaluated from the measured fraction fg of non-prompt to prompt J/y/ production in pp collisions 
and on the suppression of non-prompt J/y/ in p-Pb collisions. More in detail, in the range 

2 < ycms < 4.5, the fraction fg measured by EHCb in pp collisions at = 7 TeV, increases from 0.08 to 
0.22 from px = 0 to 8 GeV/c 1^. This quantity has a small variation within the ycms range covered and 


is also not strongly -y/ 5 -dependent (similar values are obtained for ^/s = 8 TeV |581). At mid-rapidity, /b 
was measured by ALICE in pp collisions at ^/s = 7 TeV and ranges from 0.10 to 0.44 for pj increasing 
from 1.3 to 10 GeV/c |611. R^^ prompt measured at = 5.02 TeV by EHCb, integrated over 

Pj, obtaining 0.83 ± O.Dz ± 0.08 for 2.5 < ycms < 4 and 0.98 ± 0.06 ± 0.10 for —4 < ycms < —2.5 
Assuming for each px-bin a variation of R!' 


non—prompt 
pPb 


between 0.6 and 1.3, a conservative choice due to 
the unavailability of a px-differential result, and considering the px-dependence of fg at ^/s = 7 TeV, 
one can extract as = Rp’]?)!, +/b • (^p"pb — ^ppb The maximum differences between 

the inclusive and prompt RpPb obtained in this way are, for low and high px: (i) 3 and 10% at backward 
rapidity; (ii) 11 and 16% at central rapidity; (iii) 10 and 8% at forward rapidity. These variations are, at 
most, of the same order of magnitude as the quoted uncertainties on inclusive Rppb. 


The RpPb results shown in this paper can be considered as a valuable tool to improve our understanding 
of the contribution of CNM to the suppression of the J/y/ yields observed in Pb-Pb 311. Indeed, 
as verified in |8| for the dimuon analysis, in Pb-Pb collisions the Bjorken-x ranges probed by the J/y/ 
production process in the two colliding nuclei, assuming a gg ^ J/y/ (2—)-l) mechanism, are shifted 

by only ~10% with respect to the corresponding intervals for p-Pb and Pb-p, despite the different energy 
(V'^nn = 2.76 TeV) and the slightly different ycms range (2.5 < y < 4) for Pb-Pb. A similar conclusion 
holds at mid-rapidity, where the covered v-intervals, calculated for pj = (px)> are 6.1 x 10^^ < v < 3.0 x 
10^^ and 7.0 X 10^^ < x < 3.5 x 10^^ for p-Pb and Pb-Pb collisions, respectively. Under the assumption 
that shadowing is the main CNM-related mechanism that plays a role in the J/y/ production and that its 
effect on the two colliding nuclei in Pb-Pb collisions can be factorized, the product RpPb x Rpbp (Rppb) 
can be considered as an estimate of CNM effects in Pb-Pb collisions at forward (central) rapidity | |63|[64| . 
This conclusion holds not only for the 2 —)• 1 production process but also when the more general 2 —)■ 2 
mechanism (gg — J/yfg) is considered. 


In Pig. the comparison of the measured Rpbpb with the quantities defined above is carried ouf. Such 
a comparison should be considered as qualifafive, in view of fhe slighf v-mismafch defailed above and 
of fhe facf fhaf, af mid-rapidify, fhe cenfralify ranges probed in p-Pb and Pb-Pb are nof fhe same (0- 
100% and 0-50%, respecfively). In bofh rapidify regions, fhe exfrapolafion of CNM effecfs shows a clear 
px-dependence, corresponding fo a sfrong suppression af low px, which vanishes for large fransverse 
momenfa. Af low px and cenfral rapidify, fhere mighf be an indicafion for a Pb-Pb suppression smaller 
fhan fhe CNM exfrapolafion, consisfenf wifh fhe presence of a confribufion relafed fo fhe (re)combinafion 
of cc pairs j^, faking place in fhe hof medium. A similar effecf can be seen af forward rapidify. Af 
large px and forward rapidify, fhe observed suppression in Pb-Pb collisions is much larger fhan CNM 
exfrapolafions, showing fhaf, in fhis fransverse-momenfum region, suppression effecfs in hof maffer, 
possibly relafed fo color screening, become dominant. 


Einally, a more direct comparison of Pb-Pb results with the CNM extrapolation can be obtained by 
defining fhe rafio = RpbPb/(7?pPb x ^Pbp)- Such a quantify, for forward rapidify resulfs, is shown in 


14 











Rapidity and transverse-momentum dependence of the inclusive J/yr RpPb 


ALICE Collaboration 



Fig. 6: The estimate of the px-dependence of CNM effects in Pb-Pb, calculated as /Jppj, for mid-rapidity data (top) 
and as /?ppb x /?pbp (bottom) at forward rapidity. The quantities are compared to /?pbpb measured in Pb-Pb colli¬ 
sions in the (approximately) corresponding y-ranges |30 The vertical error bars correspond to the statistical 
uncertainties, the open boxes (shaded areas) represent pp-uncorrelated (partially correlated) systematic uncertain¬ 
ties, while the boxes around Rpp\, = 1 show the size of the correlated uncertainties. The horizontal bars correspond 
to the widths of the pj bins. The Pb-Pb points in the bottom panel were slightly displaced in pj, to improve 
visibility. 
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Fig. and confirms the main features detailed above, i.e., a strong suppression of J/i/r at large pj, and 
a hint for an enhancement at low pj. At central rapidity, due to the sizeable uncertainties on both p-Pb 
and Pb-Pb results, only the px-integrated ratio can be obtained. Using the Rpbpb in the 0-90% centrality 
range |30|, and the integrated RpPb given above, one gets 1.43 ±0.26(stat) ±0.56(syst). More precise 
measurements are needed to draw a firm conclusion in this rapidity range. 



Fig. 7: The ratio between f?pbpb for inclusive J/y/ at forward rapidity and the product f?ppb x f?pbp of the nuclear 
modification factors at forward and backward rapidity. None of the uncertainties cancels out in the ratio. Statistical 
uncertainties are shown as vertical error bars, while the boxes around the points represent a quadratic combina¬ 
tion of uncorrelated and partially correlated systematic uncertainties. The box around Sj/y, — 1 corresponds to 
correlated uncertainties. The horizontal bars coincide with the widths of the pj bins. 


In summary, we have presented results on the inclusive J/i/r production in p-Pb collisions at .^/^nn = 5.02 
TeV. The px-differential cross sections, the (px) and the nuclear modification factors have been evaluated 
in three rapidity regions: —4.46 < ycms < —2.96, —1.37 < ycms < 0.43 and 2.03 < ycms < 3.53. At 
forward and mid-rapidity a significant suppression is observed at low px, with a vanishing trend at 
high Px. At backward rapidity no significant suppression or enhancement is visible. Comparisons with 
theoretical models based on a combination of nuclear shadowing and coherent energy loss effects provide 
a fair description of the observed patterns, except at forward rapidity and low transverse momentum. 
These results can be used to provide a qualitative estimate of the influence of cold nuclear matter effects 
on the J/ip suppression observed in Pb-Pb collisions. Under the assumption that shadowing represents 
the main CNM contribution, we find that it cannot account for the observed suppression in Pb-Pb at high 
Px. At low Px, the observed CNM effects alone may suggest a suppression larger than that observed in 
Pb-Pb, which is consistent with the presence of a charm quark (re)combination component to the J/ip 
production in nucleus-nucleus collisions. 
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